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ABSTRACT 


A detailed analysis of signal processing concerns for measuring aircrai 
noise is presented. Development of a de-Dopplerization scheme for both corre. 
history and spectral data is discussed along with an analysis of motion effects 
spectra. A computer code was written to implement the de-Dopplerization sch 
the code is the aircraft position data and the pressure time histories. To facilitai 
averaging, a level uniform flyover is considered in the study but the code can 
general flight profiles. The effects of spectral smearing and its removal are dis 
Using test data acquired from an XV-15 tilt-rotor flyover, comparisons are ma 
the measured and corrected spectra. Frequency shifts arc accurately accounted 
de-Dopplerization procedure. It is shown that by correcting for spherical spre; 
Doppler amplitude along with frequency can give some idea about noise sour 
The analysis indicated that smearing increases with frequency and is more sev 
approach than recession. Simulated spectra were generated using a moving pc 
model with frequency content similar to the XV-15 data. The results from the 
corroborated the observations obtained from the analysis of the actual flight d; 

INTRODUCTION 


Until recently, only 1/3 octave spectra were produced from flyover tes 
of Doppler shifts was difficult to assess with these spectra. Tonal resolution a 
frequency dependence of any broadband noise was often poor. Rapid advanc 
techniques pertaining to both algorithms and hardware have greatly influenced 
testing in recent years. This has enhanced resolution in spectra and allowed g 
flexibility in signal processing. But these advantages of digital methods requi 
knowledge of the theory behind them in order to make the most of their featu; 
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This paper describes some of the concerns and goals of signal processing 
procedures regarding aircraft flyover noise measurements. The emphasis now placed on 
HSR has produced increased interest in flyover noise measurements and its attendant data 
analysis. Since HSR developmental funding is linked to environmental issues, one of 
which is noise, accurate acoustic acquisition and signal analysis on aircraft that emulate the 
HST is of paramount importance. Due to higher speeds envisaged throughout its flight 
envelope, motion effects will be more pronounced on acoustic data for the HST than for 
conventional aircraft. 

Nonstationarity in aircraft noise and its subsequent signal restoration have been 
addressed previously. One study [ 1 ] develops two-dimensional (frequency-time) spectral 
descriptions which are based on time variant linear filter theory. This is a more accurate 
approach in theory to nonstationary analysis in the sense that it does not begin with any 
stationarity assumptions. In practice nonstationary theory increases the complexity 
involved in data analysis procedures. Nonstationary signal analysis is still in a developing 
stage and its theory is not well-known in the acoustics community. The Lorentz 
transformation is employed in another investigation |2| to relate a stationary moving source 
to its measured data which provided an estimated autocorrelation function. De- 
Dopplerization schemes that compensate for motion effects in the signal time history are the 
topics of other studies |3,4]. An approach similar to this is adopted in this paper but 
greater detail in the signal processing schemes and parameters is given along with an 
analysis of smearing effects on the measured data. 

Both the actual measured ground spectra and the de-Dopplerized versions for a 
particular aircraft are of value in determining its impact on community noise. The 
characteristics of the measured noise are what people find objectionable. But aircraft noise 
is usually quantified by static tests, thus the need for de-Dopplerized spectra. This would 
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provide some idea of a particular offending noise component where sound aba; 
should be applied. Also, the measured spectra can be used to validate current 
programs (i.e. ANOPP) and to construct exposure contours which depends o; 
microphones. In this study, a highly nonstationary signal case is considered ( 
M=.33). Using data acquired from a XV-15 Tilt-Rotor flyover test in the airpl 
assessment is made of the signal processing and de-Dopplerization techniques 

SIGNAL PROCESSING PROCEDURES 

For aircraft flyover measurements, the acoustic signal will appear nons 
ground based observer. Thus the Fourier transform of the signal will be time . 
A short-time Fourier transform with a sliding window function [5] was emplo_ 
with this. A trade-off is made here though between time resolution and freque 
resolution. One function of the window is to limit the duration of the time sigi 
spectral characteristics are reasonably stationary over the duration of the winds 
rapidly varying signal, a reduced window length is required. But reducing the 
duration also reduces the frequency resolution. Increasing the window duratic 
spectral smearing. This problem is less severe for de-Dopplerized data. 
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The detailed steps that were used in the data analysis scheme are giver; in Ref. 6. A 
brief description of these steps follows. To prevent aliasing, the bandwidth oi the 
spectrum is first selected and the sample rate determined by 


At = 


1 

2.5 f c 


( 1 ) 
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which satisfies the Nyquist criterion. Also, to avoid aliasing, the signal must be fed 
through a low-pass filter with cut-off frequency f c . Equation (1) along with the block size, 
N, determine the record length which is equal to the window duration and the bin width 
which are given by, respectively 


T = NAt 


( 2 ) 


and 



(3) 


The FFT of the signal can be expressed as |7| 


N-l 


P,(f t ) = AtP, k = At £ 

n=0 


x in w n ex H - 


i27tkn 

N 


(4) 


where Pj k is the actual output from the FFT algorithm and the subscript i designates the 
particular record. The discrete frequencies, which are the bin center frequencies, are given 
by 


f k = | , k = 0,1,2, y (5) 

With the FFT components computed, the power spectral density function can be evaluated 
from 


rid 


Gpptf^jTfll'VM 2 


i=1 


, k=0,l,2. 


N 

'2 


( 6 ) 



Using equation (6), the time-averaged pressure for the bin with center 
and width Af can be approximated by 

n d 

p2(f k ,Af) = G pp (f k )Af = ^ X |P ik (f k )| 2 


Thus, the sound-pressure level for the k-th bin is 


L(f k ) = lOlog 


p2(f k ,Af) 


Prof 


and the overall sound-pressure level for a particular band is 


Lq A = lOlog 


nb 


£ p 2 <f k .Ao 

k =o 


Prof 


MOTION EFFECTS ON MEASURED SPECTRA 

The usual model for source motion characterization is the nonhomoge 
equation i.e. 


c 2 ^ - V 2 p = q(x,t) 
° 3t2 


In this expression, the forcing function, q(x,t) , accounts for the source moth 
solution to equation (10) can be stated in terms of a volume integral [81 as 


requency f k 


(7) 


( 8 ) 


(9) 


eous wave 


( 10 ) 
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p(x,t) = — 

4k 


r 


dV(y) 


(ID 




|y-x| 


J |y 


The source is contained within the region V and y is the position vector to each source 
element. Assuming the source is concentrated at the single moving point x=x s (t) allows 
q(x,t) to be expressed as 


q(x,t) = Q(0,4>,t) 5|x-x s (t)] (12) 

where the angles 0 and <(> describe the directionality of the source. Equation (11) will then 
reduce |9] to 


p(x R ,t) 


/ r m \ Pm( 0 ><M) 

IN; 1 Mn - 


03) 


Here, p m (0,<J>,T)has been introduced to describe the source structure in terms of pressure at 
the reference distance r m . A geometrical description of this situation is shown in Figure 1. 
Equation (13) is valid for a general flight trajectory since no assumption of uniform motion 
was made in its derivation. It is important to distinguish between reception time t and 
emission or retarded time t in equation (13) which are related by 


t -x + 


M 


(14) 


In the above expressions, t is the time of signal emission and t is the time of signal 
reception. The position vector xr designates the receiver location at reception time t. Also, 
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the instantaneous Mach number vector, M, and the angles 0, <|) correspond to me x. 
Equation (14) describes the well-known signal compression and expansion di : to source 
motion. It is seen from equation (13) that amplitude modulation of the signal ill occur not 
only through r(x) and the Doppler factor (1-M n r ) but also by way of the sour e 
directionality. The next step considered is the Fourier transform of equation ( 3). 


Since p(xR,t) is the measured acoustic quantity in flyover tests, the Fo 
transform of equation ( 1 3) with respect to reception time, t, is now developed 
ascertain the effects of motion on the received spectra. The Fourier transform 
(13) is 


rier 

n order to 
»f equation 


P(x R ,f) 



p(x R ,t) e - l2n{l dt 


= r 


m 


p m (9,<l),x)e- i^ft 
|r(x)| [l-M(x)-n r (x)] 


(15) 


where the dependence on times t and x has been clearly shown since it is critic l in 
evaluating the transform. For the level flyover situation at constant velocity a shown in 
Figure 2, equation (15) becomes 


r 


P(x R ,f) = r m 


J - 


Pm(e,<M)e~ i2nfl dt 

|r(x)| (l-Mcos0) 


(16) 


Here <|) and M are constant. From the geometry, the following relation for the mission 
angle is easily derived 
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(17) 


cos0 = 


x r~ Ut 

M 


The variable of integration can be changed to the emission time x by means of equation 
(14), i.e. 


, xr-Ut 

dt = 1 1 - M-A-r-ldx 


M 


(18) 


Substituting these relations into equation (16) will result in 


P(x R ,f) = r n 


A2nf 


hi 


M 

r~ p m (0,<t>,t)e-i2nftdx 


(19) 


where r(x) can be explicitly expressed as 


r(x) = V (x R -Ux) 2 +yJ+z£ 


(20) 


Suppose the source structure is such that p m can be separated i.e. 

p m (0,<t>,t) = d(0)g(t) (21) 

where d(0) is the directivity factor and g(t) is some characteristic time signature of the 
source. Note that d can be expressed as a function of x through equation (17). Thus 
equation (19) can be written as 
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P(x R ,f) = h(t)g(t)e i2nft dx 


( 22 ) 


and h(x) is defined as 


- ErfN 

h(x)=% rM(x) 

m\ 


(23) 


Therefore, P(xr,0 is the convolution of the Fourier transforms of h(t) and g(t , namely 


t ,f)=j 


P(x R ,f) = H(a)G(f-a)da 


(24) 


As a final consideration of some of the theoretical aspects of ground-fc sed spectra, 
approximate relations for a short-time Fourier transform will be given. The at ;uracy of 
these approximations is very much dependent on source velocity. Assuming ■ tat t is short 
enough such that Ux « IrRl permits the following approximations to be made 


M * f R 


Ux R t 

tr 


(25) 


and 


cos 0 = 


XR 

r R 


(26) 


Substituting these expressions into equation (19) results in 
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i r R f 
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P(x R ,f) = -yV p m (6,<)>,x)e i2nr ( i-Mcose)T d x 


(27) 


If the smear angle, A9, is small, a characteristic source spectrum can be defined as 


( 


P m (0,<t>,f) = I p m (0,<j),x)e i 2 nft dx 


(28) 


So that equation (27) can be expressed as 


,2n( 

-i r K 


P(x R ,f) = P m |0,<>,f(l-Mcos0)] 

r R 


(29) 


For a pure tone of amplitude A and frequency f 0 , Pm is 


sin[7t(f-f 0 )T] 

P m (0,<J>,f) = A T e- i*r-UT 

7C(f-f 0 )T 


(30) 


which is the well-known result describing spectral broadening due to a finite time duration. 
Introduction of equation (30) into equation (29) produces 


n/ ^ A^e-iP sin{n|f(l-Mcos0)-f o ]T) 

P(x R ,0 = AT — 

r R n[f(l-Mcos0)-f o lT 


(31) 


The phase angle (3 is given by 



( 32 ) 


(3 = ^r R + jt[f(l -Mcos0)-f o ]T 
C n 


From equation (31), it is seen that its first zero occurs at 


f, = 


1 r 
^ 

1-McosG 


(33) 


For an approaching source (subsonic case) 1-Mcos0<l and for a recec ng source 1- 
Mcos6>l. Thus, fi is greater in the approach regime and this indicates spectr; smearing is 
more pronounced in this portion of the flight path. 


DE-DOPPLERIZATION SCHEME 


(34) 

Taking the Fourier transform of equation (34) with respect to x will result in tl ; source 
spectrum with corrections for amplitude and frequency, viz. 


Equation (13) can be re-written as 


p m (0,<t>,x) = 


Ml 


l m / 


(l-M n r )p(x R ,t) 


Pm(0’<M) = I - M-n r )p(x R ,t)e '2*fTdi 

I * m 


(35) 


So that, in principle p m can be evaluated using the measured data. Th< tracking 
data determines r(x) and M while the measured acoustic data provides p(xR,t) Equation 
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(14) allows p to be expressed as a function of the emission time x. The specific steps used 
in this study to compute equation (35) were as follows. First, at a particular microphone 
location a start time is chosen to begin the analysis. Then the tracking data is searched for 
this time where linear interpolation is used to determine r(x) and M(x). That is, this start 
time is emission time. The acoustic data, p(xR,t) , is in a digital format consisting of 
equally-spaced samples. Equation (14) is employed to compute the reception time t that 
corresponds to the emission time x. Linear interpolation is performed on the p vs. t time 
history to determine p(xR,t). Incrementing x by the designated sample rate At generates an 
equally-spaced pressure time history which is the required form for the FFT algorithm. 
These steps are summarized in Figure 3. A word of caution concerning aliasing in this 
signal reconstruction process. For an approaching aircraft, the measured signal is 
compressed with the highest frequency limited to f c by the anti-aliasing filter. The 
reconstructed signal will appear to be expanded and its highest frequency will be less than 
f c . Thus, there is no problem with aliasing. But for a receding aircraft the measured signal 
is expanded and the restored signal will appear to be compressed. Therefore it could have 
spectral content above f c . In this situation aliasing could occur. If significant power is 
expected above f c then the reconstructed signal should be low-pass filtered or re-sampled at 
a higher rate. 

So far in the above discussion no stipulation has been made for a level, constant 
velocity flyover. To implement an ensemble average the uniform, level flyover as shown 
in Figure 4 is considered. An ensemble average is performed across four equally-spaced 
measurement locations beneath the flight path. The records used in the average are taken at 
approximately the same emission angles and smear angles. A description of the test and the 
data acquisition for the XV- 15 is given in Ref. 10. 
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Rush mounted microphones were employed which had a frequency re 
0.01-40,000 Hz with a maximum sound-pressure level of 160 dB. The signal 
passed filtered between 16-16,000 Hz and then RVI -recorded at 15 ips which j 
dynamic range of about 46 dB. The recorded signals were then low-pass filter 
kHz and digitized at a sample rate of 25 kHz. All of the results presented in th 
developed from one flyover measurement of the XV- 15 in the airplane mode ( 
M=.33). The aircraft track is shown in Figure 5 which is determined by a lase 
system. Since the sample rate is set, the only way to vary the bin width is thro 
block size N. For the first group of results discussed N is set to N=2048. Th 
corresponds to a window duration of T=.08192s and a bin width of Af=12.2 i 
averaging was performed within each microphone measurement and then the n 
spectra averaged across the four microphones 1 10]. At each measurement loci 
acquired time history was divided into five contiguous records. Thus in these 
Spectra were produced from seven measurement intervals of the flight path. T 
contains the averaged emission angles, smear angles and frequency shifts for t 
The frequency parameters, fR and AfR are calculated from the following 
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(36) 


and 


Af„ = 




1 1 _ 

1-Mcos0| 1-Mcos02 


(37) 


which are the Doppler shift relations for uniform motion. Equation (37) gives i 
quantitative idea about spectral smearing for a window duration that spans the mear angle 

given by 
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A0 = 0 2 - 0i 


(38) 


Figures 6 and 7 show the time histories that correspond to microphone no.l for 
©1=16.25°. In Figure 6, the measured pressure is illustrated. Figure 7 depicts the 
corrected signal using equations (14) and (34). Here, r m =62.5 ft. which is 5x the propeller 
radius of the XV- 15. Also, the plot in Figure 7 was constructed using equally-spaced 
samples. The signal compression in the measured data is clearly seen by comparing the 
figures. For this test flight, the propellers were operated at 589 rpm which yields a blade 
passage frequency of 29.45 Hz. Notice that an effective period of .039s can be ascribed to 
the signal in Figure 7 which yields a fundamental of 25.64 Hz. Measurements to the rear 
of the aircraft are presented in Figures 8 and 9 which were taken by microphone no.l for 
01=141.15°. These are similar to the previous figures in that Figure 8 is the measured 
signal and Figure 9 is the de-Dopplerized signal. Now, expansion of the signal in the 
measured data is evident by comparing the figures. Though the time histories differ in 
structure between fore and aft measurement locations they do exhibit the same characteristic 
period of ,039s. Thus the de-Dopplerization procedure is accurately accounting for signal 
compression and expansion. 

Figures 10-16 describe the measured (shifted) spectra for a bandwidth of 10 kHz. 
These spectra confirm 1 10) that most of the power and structure of the signal in the airplane 
mode is concentrated below 1 kHz. A dominant tone occurs in the neighborhood of the 
blade passage frequency. The de-Dopplerized spectra are shown in Figures 17-23. The 
apparent tonal structures at 8.56 kHz in Figure 17, 8.7 kHz in Figure 18 and 9 kHz in 
Figure 19 are artifacts due to the low-pass filter. In each of these cases, these frequencies 
yield a Doppler shifted frequency of 12.5 kHz, which is the Nyquist frequency for the 
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data. Thus, for these figures, the portion of the spectra above these respet ive frequencies 
should be ignored since they represent the transition zone of the filter. 

The spectra in Figures 17-20 represent radiation to the fore of the p opeller plane. 
They were determined from a sweep of 0 through the range 16.25°-65.03° The four 

spectra show little change. There is a slight increase in the fundamental as he propeller 
plane is approached. Table 2 contains the peak and overall sound-pressure levels for all 
seven blocks for both Doppler and de-Dopplerized spectra. In the fore reg on of the 
propeller, it appears that the directivity pattern remains relatively uniform tl rough this range 
of 0. Figures 22 and 23 depict sound propagation to the aft of the propelle plane. These 

spectra and Table 2 indicate greater variability in the radiation pattern in the aft region 
compared to the fore region. Since the peak and overall sound-pressure le 1 els are less in 
this portion of the acoustic field, suggests that the energy radiated to the aft of the aircraft is 
less than to the front. Also, as is the case in the fore direction, the levels d< ;rease as the 
propeller axis is approached. 

After this analysis an attempt to improve the resolution of the specti i by increasing 
the block size, N, was performed. Choosing N= 16384 produces a windov duration of 
T=.65536s and a bin width of Af=1.53 Hz. Due to the increased window . uration no 
averaging was attempted within each microphone measurement but only at oss the four 
acquisition locations so that n ( p4. As before, spectra were constructed froi i seven 
measurement intervals of the flight path. Table 3 shows the greater smear ; lgles that result 
for this increased resolution. Figures 24-30 contain the spectra of the meas ired signals for 
these cases for a bandwidth of 500 Hz. Likewise, Figures 31-37 contain th : de- 
Dopplerized spectra. Table 4 presents the sound-pressure levels. The fact hat for some 
cases in Table 4 the overall levels are greater than those shown in Table 2 e en though the 
frequency range is less can be attributed to the finer frequency resolution, i l the Doppler 
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spectra smearing is evident beginning with the 4th harmonic in Figure 26 and the 2nd 
harmonic in Figures 27-29. Also note the broadening of the fundamental in these figures. 

In general, the first seven harmonics are easily identified in the de-Dopplerized 
spectra. Figures 31-37. The pronounced smearing that was evident in the first three 
overtones in the measured spectra is absent in these figures. Though the amplitudes of the 
harmonics vary during the flyover, they do retain their spectral location in each spectrum. 
The fundamental occurs at approximately 25 Hz. As expected, since most of the energy is 
contained in the lower end of the spectrum, the same trends emerge that were produced for 
a 10 kHz bandwidth and Af= 12.21 Hz. That is, the levels increase toward the propeller 
plane with more energy radiated to the front of the aircraft and greater variation in the 
radiation pattern to the rear of the aircraft. An important feature brought out by the finer 
resolution of this data is the difference in the fundamental and the overtones in the fore and 
aft directions. As Figures 31-34 show the fundamental dominates the peak overtone by 18- 
21 dB in the forward direction. But, Figures 36 and 37 show that to the rear of the aircraft 
the 2nd harmonic is approximately only 2 dB less than the fundamental. 

SIMULATED SPECTRA 

To complete the study, an analysis of some computer simulated spectra is now 
presented. With known source distribution and flight path, a pressure time history can be 
constructed at the observation point by way of equation (13). Therefore, a simulated 
spectrum is generated at the measurement location of the moving source. All the cases 
depicted were produced using a level flyover (h=250 ft.) and U=222.3 knots (M=.33). 

The idea being to mimic as closely as possible the situation in the previous section. An 
omnidirectonal source was chosen whose frequency content consisted of a summation of 
harmonics. Thus the form of p m (0,<J>,x) is given by 
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N„ 

Pm( 0 ><M) = X A n co 
n=l 

where 

f n = nf i ,n=l,2 N h (40) 

This yields a periodic signal with a fundamental frequency of fj. I y choosing 
f 1 =25 Hz and Aj equivalent to a 115 dB tone (re 62.5 ft.) yields similar ch racteristics to 
the real data. The amplitudes of the overtones monotonically decreased wi i each harmonic 
25 percent less than the previous one (A n =.75A n .|) and the number of han ionics in the 
signal was set to Nh=100. As before time series were generated in seven I locks of the 
flight path where each block began at the same emission angles given in T bles 1 and 3. 
Since the frequency content of the source is known it is easily seen how w< il the signal 
analysis scheme is performing. 

Using altitude, velocity and initial emission angle, 0o, as input, the x>sition vectors 
can be expressed as 

x s Cc) = x s i+ysj+z s k = Ut i (41 ) 

and 

X R = x R i+ yiJ +z R k = hcotSoi - hj (42) 

Also, equation (14) can be solved explicitly for x in terms of t 


2*f„t+{n-l)j 
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( 43 ) 


c 0 t - Mx r - V(x R -UtJ2+(l-M2)y2 
X " c 0 (l-M2) 

Equation (43) is needed for signal processing based on reception time (contains 
Doppler shifts) since FFT algorithms require equally-spaced samples. For the results now 
presented, a bandwidth of 500 Hz and a bin width of Af=1.53 Hz was again chosen. 
Figure 38 illustrates the pressure time history based on emission time which is given by 
equation (39). This can be compared to the received signals computed from equation (13) 
shown in Figures 39 and 40. Signal compression for an approaching source is illustrated 
in Figure 39 for 01=16.25°. Signal expansion for a receding source is depicted in Figure 
40 for 01=141.15°. Figures 41-47 contain the spectra of the seven data blocks and Table 5 
shows some of the pertinent characteristics of the spectra. The increase in spectral 
smearing with frequency is evident. It is more severe on approach becoming apparent on 
the lower harmonics as the smear angle increases. From Figure 41 it is seen that tonal 
broadening first becomes significant at the 7th harmonic. But Figure 47 which was 
constructed as the source receded (01=141.15°) shows that the harmonics are discernible 
through the 1 1 th harmonic. On approach, notice how the smearing of the higher-order 
harmonics takes on the appearance of band-pass white noise. 

CONCLUSIONS 

A detailed analysis of signal processing concerns for measuring aircraft flyover 
noise has been presented in this paper. The effects of spectral smearing and its removal 
(de-Dopplerization) were considered. Using test data acquired from an XV-15 tilt-rotor 
flyover, comparisons were made between the measured and corrected spectra. Frequency 
shifts are accurately accounted for by the de-Dopplerization scheme. It was shown in the 
study that by correcting for spherical spreading and Doppler amplitude along with 
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frequency can give some idea about noise source directivity. Also, the analys ; indicated 
that smearing increases with frequency and is more severe on approach than i cession. 


Simulated spectra were generated using a moving point source model with fre uency 
content similar to the XV- 15 data. The results from the simulation corroborat d the 
observations obtained from the analysis of the actual flight data. 
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APPENDIX: PROGRAM LISTING 


PROGRAM BBAND 

IMPLICIT REAL*8 (A-H,0-Z) 

REAL STI(32766) , SPR(40000) , SFREQ(2048) , SSR(2048) 

REAL SIR(50000) , SIT(50000) 

DIMENSION PR( 32766 ) , SR(2048) , FREQ (2048) , IWK(20) 

DIMENSION TI(32766) , SRI (2048) , TSO(4 , 9) 

DIMENSION TR( 50000) , PRE(50000) , PRT(50000) 

DIMENSION TIM(2048),X(2048),Y(2048),Z(2048) 

DIMENSION VX( 2048 ) , VY( 2048 ) , VZ( 2048 ) , NRD( 4 ) 

DIMENSION SX(2048),SR1X(2048) 

COMPLEX*16 ZS( 16384), ZX( 16384) 

OPEN(UNIT=20 , FILE= ' FLYOR . OUT ' , STATUS= ' NEW' ) 

OPEN ( 8 , FI LE= ' KRAD0200 . 909 ' , F0RM= ' FORMATTED ' , ACCESS= 

+ ' SEQUENTIAL' , STATUS=' OLD' ) 

************************************************************* 
OPEN ACOUSTIC FILES 

************************************************************* 

IBLK=24000 

IBUF=1 

TS0(1, 1 ) =36349. 7006D0 
TSO(2 , 1 ) =36350. 2354D0 
TS0(3 , 1 ) =36350. 77 16D0 
TSO(4, 1) =36351 . 3041D0 
OPEN ( 1 , FI LE= ' KA02009A . DAT ' , 

+ FORM= ' UNFORMATTED ' , 

+ RECORDTYPE=' FIXED' ,BLOCKSIZE=IBLK,BUFFERCOUNT=IBUF, 

+ ACCESS=' DIRECT' ,RECL=2048,ASS0CIATEVARIABLE=IPLACE, 

+ STATUS='OLD' ) 

0PEN(2 , FILE= ' A020010A.DAT' , 

+ FORM= ' UNFORMATTED ' , 

+ RECORDTYPE=' FIXED' ,BLOCKSIZE=IBLK,BUFFERCOUNT=IBUF, 

+ ACCESS= ' DIRECT ' , RECL=2048 , ASSOCI ATEVARIABLE=IPLACE , 

+ STATUS=' OLD' ) 

0PEN(3, FILE=' A020011A.DAT' , 

+ F0RM= ' UNFORMATTED ' , 

+ RECORDTYPE=' FIXED' , BLOCKSIZE=IBLK, BUFFERCOUNT=IBUF, 

+ ACCESS=' DIRECT' , RECL=2048 , ASSOCIATEVARIABLE=IPLACE , 

+ STATUS='OLD' ) 

0PEN(4,FILE=' A020012A.DAT' , 

+ FORM= ' UNFORMATTED ' , 

+ RECORDTYPE=' FIXED' , BLOCKSIZE=IBLK, BUFFERCOUNT=IBUF, 

+ ACCESS=' DIRECT' , RECL=2048 , ASSOCIATEVARIABLE=IPLACE , 

+ STATUS=' OLD' ) 

39 F0RMAT(/,1X, 'NREC=' ,15) 

DO 300 K=l,4 

CALL HEADER_BUFFER ( K , NREC ) 

WRITE(20, 39)NREC 
300 NRD(K)=NREC 

READ RADAR DATA 

TSTART=1ST ACOUSTIC TIME MARK 

TSTART=36349 . 7006D0 
MS=1 

1 READ(8,101,END=654)TIM(MS) ,X(MS),Y(MS),Z(MS), 

+VX(MS) ,VY(MS) ,VZ(MS) 

101 F0RMAT(1X, Fll . 1 , 6F10. 1 ) 

IF(TIM(MS) . LT. (TSTART- . 3D0) )G0 TO 1 
X(MS)=-X(MS) 

VX(MS)=-VX(MS) 21 

MS=MS+1 



non 


GO TO 1 
654 MS=MS-1 

0 ************************************************** 

C SIGNAL PROCESSING PARAMETERS 

C FC=BANDWIDTH, NPTS=N0. OF SAMPLES, NR=NO. OF CHANNELS 
C DT=SAMPLE RATE, TO=TEM PERATURE ( C ) 

C XR, YR, ZR=MICROPHONE LOCATION 
0 ************************************************** 
FC=500.D0 
M=14 

NPTS=2**M 
NR=4 

COEF=2 . DO/DFLOAT(NR*NPTS**2 ) ! BOXCAR WINDOW 
DT=4.E-5 

PREF=2 . E-4 ! SIGNAL IN DYNES/CM**2 
RREF=62 . 5DO !REF. DISTANCE IN FT 
IP=0 
IX=0 
IY=0 

BW=1 . DO/ ( DT*DFLOAT(NPTS) ) 

TO=21 . 9DO 
YR=0 . DO 
ZR=0 . DO 

WRITE (20, 50 )TO 

50 F0RMAT(//,1X,'T0=' ,F5.1,'C' ,//) 

C0=DSQRT(401 .8DO*(TO+273. 16D0) )/ . 3048D0 
WRITE(20,60)DT,BW,NR,MS 

60 FORMAT(// , IX, ' DT=' , E9 . 3 , ' SEC' , 2X, ' BW=' , F5. 2 , ' HZ' , 2X, 
* ' NR= ' ,13, 2X, ' #PTS=' ,15,//) 

LI=FC/BW 
DO 25 Kl=l , LI 
FREQ(K1)=DFL0AT(K1-1)*BW 
25 SFREQ(K1)=SNGL(FREQ(K1>) 

CALL PSEUDO 
CALL PSE(ll) 

GENERATE DIGITAL INPUT SIGNAL 

DO 5 NT=1 , 8 
DO 20 LM=1 , LI 
SX(LM)=0. DO 
20 SR(LM)=0. DO 
SUM=0 . DO 
SUMX=0 . DO 
XR=0 . DO 
NSS=32765 

IF(NT.GE. 5)NSS=32765 
IF(NT.EQ.8)NSS=22730 
IMEML0C=1 + (NT- 1) *12500 
DO 6 ICHAN=1,4 

CALL DATA_BUFFER( ICHAN , NRD( ICHAN) , IMEMLOC , NSS , SPR ) 
TSO(ICHAN,NT)=TSO(ICHAN, 1 )+. 5D0*DFL0AT(NT-1) 

DO 45 NN=1 ,NSS 
PR(NN)=DBLE(SPR(NN) ) 

45 TI(NN)=TSO( ICHAN, NT) +DFL0AT(NN-1)*DT 
NX=1 
NY=1 

DO 51 MX=1 ,NPTS 
IF(TI(MX) .GE.TIM(MS) )STOP 

CALL NTRP0L8VEC( X , Y , Z , VX , VY , VZ , TIM , MS , TI ( MX) , 

+RX , RY , RZ , RVX , RVY , RVZ , NX , MS ) 

RDX=XR-RX 
RDY=YR-RY 
RDZ=ZR-RZ 

R=DSQRT(RDX**2+RDY**2+RDZ**2 ) 22 

UX=RDX/R 
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UY=RDY/R 

UZ=RDZ/R 

DP=UX*RVX+UY*RVY+UZ*RVZ 

DEN0M=1. DO-DP/CO 

IF(MX. EQ. 1 . OR. MX. EQ.NPTS)THEN 

THE=DAC0SD(DP/DSQRT(RVX**2+RVY**2+RVZ**2) ) 

WRITE (20, 75 )THE, I CHAN 
75 FORMAT (/ , IX , ' THETA= ' , F5 . 1 , 2X , ' MIC= ' , 12 ) 

ENDIF 

TR(MX)=TI(MX)+R/C0 

PRT(MX)=0.D0 

PRE(MX)=0.D0 

IF(TR(MX) .GE.TI(NSS) )G0 TO 51 

CALL NTRP0L8(PR,TI, 32766, TR(MX), RESULT, NY, NSS) 

PRT(MX)=RESULT 

PRE( MX ) =RESULT* ( R/RREF ) *DEN0M 
51 CONTINUE 

PLOT TIME HISTORY (MEASURED SIGNAL) 

IF(IP.EQ.0)G0 TO 4 
DO 100 11=1, NSS 
100 STI(II)=SNGL(TI(II )-TI( 1) ) 

CALL PL0TR( 1 , 0 , NSS , STI , SPR , 0 . ,0. ,0. ,0. , 'TIME(SEC) ' ,9, . 16 , 
*'P(DYNES/CM**2)' ,14, .16, 

*' MEASURED TIME HISTORY' , 21 , . 16 ,0,0. , 10. , 1 . 5 , 13. 6 , 7 . , 1 . , 
*10.4,0,1) 

CALL ERASE 

PLOT TIME HISTORY (RECEPTION TIME) 

IF(IX.EQ.0)G0 TO 4 

DO 200 11=1, NPTS 

SIT(II )=SNGL(TR(II)-TR( 1 ) ) 

200 SIR( II )=SNGL( PRT ( II ) ) 

CALL PLOTR( 1,0, NPTS, SIT, SIR, 0. ,0. ,0. ,0. , 'TIME(SEC)' ,9, . 16, 
*'P(DYNES/CM**2)' ,14, .16, 

*' RECEPTION TIME HISTORY' , 22 , . 16 , 0, 0 . , 10. , 1 . 5, 13. 6, 7 . , 1 . , 
*10.4,0,1) 

CALL ERASE 

PLOT TIME HISTORY (EMISSION TIME) 

IF(IY.EQ.O)GO TO 4 
DO 205 11=1, NPTS 
205 SIR(II)=SNGL(PRE(II) ) 

CALL PLOTR( 1,0, NPTS, STI, SIR, 0. ,0. ,0. ,0. , 'TIME(SEC) ' ,9, . 16, 
*'P(DYNES/CM**2)' ,14, .16, 

★'EMISSION TIME HISTORY' , 21 ,. 16,0,0. , 10. , 1 . 5, 13.6,7 ., 1 . , 
*10.4,0,1) 

CALL ERASE 

COMPUTE SPL FOR EACH BIN AND OASPL 

4 DO 2 J=1,NPTS 

ZX(J)=DCMPLX(PR( J) ,0.D0) 

ZS( J)=DCMPLX(PRE( J) ,0. DO) 

2 CONTINUE 

CALL FFTD(ZX,M, IWK) 

CALL FFTD(ZS,M, IWK) 

DO 3 JK=1,LI 
SRX=(CDABS(ZX( JK) ) )**2 
SRT=(CDABS(ZS( JK) ) )**2 
SX( JK)=SX( JK)+SRX 

3 SR( JK)=SR( JK)+SRT 23 

6 XR=XR-200. DO 



DO 30 LS=1 , LI 
SR1X(LS)=C0EF*SX(LS) 

SUMX=SUMX+SR1X( LS) 

SR1(LS)=C0EF*SR(LS) 

SUM=SUM+SR1(LS) 

SX(LS)=10. D0*DL0G10(SR1X(LS)/PREF**2) 

30 SR( LS ) =10 . D0*DL0G10( SRI ( LS) /PREF**2 ) 

OASPL=10 . D0*DL0G10( SUM/ PREF**2 ) 

OASPLX= 10 . D0*DL0G10(SUMX/PREF**2) 

C 

C PLOT SPECTRUM (DOPPLER SPECTRUM) 

C 

DO 105 11=1 , LI 

105 SSR(II)=SNGL(SX(II) ) 

CALL PL0TR( 1 , 0, LI , SFREQ, SSR, 0. ,0. ,0. ,0. , ' FREQ(HZ) ' ,8,-16 
* ' SPL( dB ) ' , 7 , . 16 , ' DOPPLER SPECTRUM ' , 16 , 

*.16,0,0. ,10. ,1.5,13.6,7. ,1. ,10.4,0,1) 

CALL ERASE 
C 

C PLOT SPECTRUM (DE-DOPPLERIZED SPECTRUM) 

C 

DO 106 11=1, LI 

106 SSR(II )=SNGL( SR(II ) ) 

CALL PLOTR( 1,0, LI, SFREQ, SSR, 0. ,0. ,0. ,0. , 'FREQ(HZ)' ,8, .16 
*' SPL(dB) ' , 7 , . 16 , ' DE-DOPPLERIZED SPECTRUM' , 23 , 

*.16,0,0. ,10. ,1.5,13.6,7. ,1. ,10.4,0,1) 

CALL ERASE 

WRITE( 20 , 70)0ASPLX, OASPL, NT 

70 FORMAT(/ , IX, ' OASPLX=' , F4.0, ' dB' ,2X, 'OASPL=' ,F4.0, 'dB' ,2X 
*'NT=' ,12) 

5 CONTINUE 

CALL CALPLT(0. ,0. ,999) 

CALL PSC 
CLOSE(20) 

STOP 

END 

C 

C INTERPOLATION SUBROUTINE 
C 

SUBROUTINE NTRP0L8( VC1 , VC2 , N , REF , RESULT , I , M ) 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION VCl(N) , VC2(N) 

IKEEP=0 

1 IF( (VC2(I )/REF) . LE. 1 . DO)THEN 
IKEEP=I 

I F ( I KEEP . EQ . M ) STOP 
1=1+1 
GO TO 1 
ENDIF 
1=1-1 

IF( IKEEP . EQ . 0 ) STOP 
V2DIFF=VC2 ( IKEEP+1 ) -VC2 ( IKEEP) 

V1DIFF=VC1( IKEEP+ 1 ) -VC1 ( IKEEP ) 

RESULT=(V1DIFF/V2DIFF)*(REF-VC2 ( IKEEP+1 ) )+VCl (IKEEP+1) 
RETURN 
END 
C 

C RADAR DATA INTERPOLATION 
C 

SUBROUTINE NTRP0L8VEC ( VC1 , VC2 , VC3 , VC4 , VC5 , VC6 , VC7 , 

+N , REF , RES 1 , RES2 , RES3 , RES4 , RES5 , RES6 , I , M ) 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION VC1(N) , VC2(N) ,VC3(N) , VC4(N) , VC5(N) 

DIMENSION VC6(N) , VC7(N) 

1 IF((VC7(I)/REF).LE.1.D0)THEN 24 

IKEEP=I 



IF(IKEEP.EG.M)STOP 
1=1 + 1 
GO TO 1 
ENDIF 
1 = 1-1 

IF( IKEEP. EQ. 0)ST0P 
V7DIFF=VC7 ( IKEEP+1 ) -VC7 ( IKEEP) 
V6DI FF=VC6 ( IKEEP+ 1 ) -VC6 ( IKEEP ) 
V5DI FF=VC5 ( IKEEP+ 1 ) -VC5 ( IKEEP ) 
V4DIFF=VC4( IKEEP+1) -VC4 (IKEEP ) 
V3DIFF=VC3(IKEEP+1 )-VC3( IKEEP) 
V2DIFF=VC2 ( IKEEP+1 ) -VC2 ( IKEEP ) 
V1DIFF=VC1 (IKEEP+1 )-VCl (IKEEP) 
DTA= ( REF-VC7 ( IKEEP+ 1 ) ) / V7DIFF 
RES1=V1DIFF*DTA+VC1(IKEEP+1 ) 
RES2=V2DIFF*DTA+VC2 ( IKEEP+1 ) 
RES3=V3DIFF*DTA+VC3 ( IKEEP+ 1 ) 
RES4=V4DIFF*DTA+VC4( IKEEP+1) 
RES5=V5DIFF*DTA+VC5( IKEEP+ 1 ) 
RES6=V6DIFF*DTA+VC6(IKEEP+1 ) 
RETURN 
END 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

*= 


*- 


SUBROUTINE TO EMULATE CALLING MECHANISM TO GRANDLE'S A/D BUF1 2RS 


VARIABLE 

RANGE 

DESCRIPTION 

I CHANNEL 

0-5 

(WHICH MIC-CHANNEL TO USE) 

THERE WILL CONCEPTUALLY BE 
FIVE ROUTINES LIKE THIS , EACH 
CAN ACCESS ONLY 6 CHANNELS. 

THIS IS THE CASE SO AS TO MIMMICK 
THE HARDWARE CONSTRAINT. 

IFIRST 

0,16777215 

INTEGER ARGUMENT INDICATING WHICH 
ELEMENT IN THE DATA STREAM TO BE 
ACCESSED. 

IHOWMANY 

0,32765 

AT ANY ONE CALL THIS SUBROUTINE 
WILL RETURN THIS MANY ELEMENTS 
IN THE DATA ARRAY. 

DATA 


ARRAY WHICH CONTAINS DATA. 


SUBROUTINE header buffer ( i channel , nor i f ) 

j header read ei try point 

READ(ICHANNEL'l) NORIF ! NUMBER OF RECi RDS IN FILE 

return ! return to mail prog 

end 


subroutine da t abuf fer ( i channel , nrec , i first , ihowmany , data) 

* = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = __ = = = = = = = = = = _f ( j ata rea( J ent , y p 0 i nt 

dimension data(40000) 
data nvpr / 2048 / 


IFREC2RD = INT(IFIRST/NVPR) + 2 IFIRST RECORD I) READ FROM 

!THE FILE 
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ILREC2RD = INT((IFIRST~1+IH0WMANY)/NVPR) + 2 !LAST RECORD TO READ FROM 

!THE FILE 

NUMRECS = ILREC2RD + 1 - IFREC2RD ! TOTAL NUMBER OF RECORDS 

!TO READ 


* issue error diagnostic if attempt to read beyond file limits 

★ 

★ 

IF( ( IFREC2RD + NUMRECS - 1) .GT. nrec) THEN 

TYPE ' ATTEMPT TO READ BEYOND EOF CHANNEL READ FAILED' 
RETURN 

ENDIF 


* READ THE PERTINENT DATA FROM THE DIRECT ACCESS FILE INTO THE DATA ARRAY 


DO irecp = IFREC2RD, ILREC2RD 
L = K + NVPR - 1 

READ (I CHANNEL' irecp) (DATA( J) , J=K,L) 
K = K + NVPR 

ENDDO 


★ 

* THROW AWAY DATA IN THE BEGINNING OF THE ARRAY THAT WAS NOT REQUESTED 

★ 

* 


IFRSTNREC = (IFREC2RD - 2)*NVPR 10-2047 INDEX SCHEME ID FOR THE FIRST 

! ELEMENT OF THE DATA ARRAY 

INC = IFIRST - IFRSTNREC - 1 ! SHIFTING INCREMENT CONSTANT 

* 

* fill data array with requested data 

* 

★ 

DO I = 1 , IHOWMANY 
DATA(I) = DATA (INC + I) 

ENDDO 


RETURN 

END 
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Table 1 

Emission angles and Doppler shifts for 
Af =12.2 Hz 


block no. 

0 1 (deg) 

e 2 (deg) 

A0 (deg) 

Af R 

1 

16.25 

19.25 

3.0 

.001 

2 

20.3 

25.98 

5.68 

.026 

3 

27.6 

38.53 

10.93 

.065 

4 

41.9 

65.03 

23.13 

.164 

5 

72.0 

107.78 

35.78 

.205 

6 

114.7 

137.73 

23.03 

.075 

7 

i 

141.15 

152.1 

10.95 

.021 
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Table 2 



Doppler spectrum 

de-Dopplerized spectrum 

peak spl 

m 

oaspl (dB) 

peak spl 
(dB) 

oaspl (dB) 

90 

91 

114 

115 

94 

95 

115 

116 

98 

99 

116 

117 

103 

104 

117 

118 

102 

107 

115 

117 

100 

102 

109 

113 

90 

94 

103 

108 
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Table 3 

Emission angles and Doppler shifts for 
Af =1.53 Hz 


block no. 

0 1 (deg) 

e 2 (deg) 

A0 (deg) 

1 

< 

1 

16.25 

22.5 

6.25 

.025 

2 

20.3 

31.03 

10.73 

.054 

3 

27.6 

48.93 

21.33 

.13? 

4 

41.9 

85.25 

43.35 

.16-' 

5 

72.0 

125.1 

53.1 

1 

^ j 

o 

CNJ 

6 

114.7 

146.1 

31.4 

.09- 

7 

141.15 

156.25 

15.1 

.02 > 
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Table 4 


sound-pressure levels 
for Af =1.53 Hz 
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Table 5 

Emission angles and sound-pressure levels for 
simulated spectra, Af = 1 .53 Hz 


block no. 


0i (deg) 

0 2 (deg) 

A0 (deg) 

peak spl 
(dB) 

o ispl (dB) 

16.25 

26.51 

10.26 

96 

101 

20.3 

37.6 

17.06 

98 

103 

27.6 

59.05 

31.45 

98 

105 

41.9 

91.95 

50.05 

97 

107 

72.0 

121.69 

49.69 

97 

106 

114.7 

141.82 

27.12 

95 

102 

141.15 

153.54 

12.39 

95 

99 
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Figure 1 : Flight geometry and parameter description 
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Figure 2: Level flyover geometry 



Figure 3. Time history restoration scheme 












Figure 5: Nasa/Bell XV-15 position data - level flyover run 




Figure 6. Measured time history; microphone no. 1; 01 = 16.25°; 
A0=3.O°. 
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Figure 7. Emission time history; microphone no. 1; 01 = 16.25°; 
A6=3.0°. 
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Figure 11. Doppl 
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Figure 22. De-Dopplerized spectrum; Af=12.2 Hz; 0 i = ll 4 . 7 O ; 
A0=23.O3°. 
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42. Simulated spectrum; Af=1.53 Hz; 01=20.3°; A0-17.O6 0 . 
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